A mathematical model of ENaC and Slc26a6 regulation by CFTR in salivary gland ducts.

Cystic fibrosis (CF) is a genetic disease caused by the mutations of cystic fibrosis transmembrane conductance regulator (CFTR), the cystic fibrosis transmembrane conductance regulator gene. Cftr is a critical ion channel expressed in the apical membrane of mouse salivary gland striated duct cells. Although Cftr is primarily a Cl- channel, its knockout leads to higher salivary Cl- and Na+ concentrations and lower pH. Mouse experiments show that the activation of Cftr upregulates epithelial Na+ channel (ENaC) protein expression level and Slc26a6 (a 1Cl-:2[Formula: see text] exchanger of the solute carrier family) activity. Experimentally, it is difficult to predict how much the coregulation effects of CFTR contribute to the abnormal Na+, Cl-, and [Formula: see text] concentrations and pH in CF saliva. To address this question, we construct a wild-type mouse salivary gland model and simulate CFTR knockout by altering the expression levels of CFTR, ENaC, and Slc26a6. By reproducing the in vivo and ex vivo final saliva measurements from wild-type and CFTR knockout animals, we obtain computational evidence that ENaC and Slc26a6 activities are downregulated in CFTR knockout in salivary glands.NEW & NOTEWORTHY This paper describes a salivary gland mathematical model simulating the ion exchange between saliva and the salivary gland duct epithelium. The novelty lies in the implementation of CFTR regulating ENaC and Slc26a6 in a CFTR knockout gland. By reproducing the experimental saliva measurements in wild-type and CFTR knockout glands, the model shows that CFTR regulates ENaC and Slc26a6 anion exchanger in salivary glands. The method could be used to understand the various cystic fibrosis phenotypes.

Tools for Quantitative Analysis of Calcium Signaling Data Using Jupyter-Lab Notebooks.

Calcium signaling data analysis has become increasing complex as the size of acquired datasets increases. In this paper we present a Ca2+ signaling data analysis method that employs custom written software scripts deployed in a collection of Jupyter-Lab "notebooks" which were designed to cope with this complexity. The notebook contents are organized to optimize data analysis workflow and efficiency. The method is demonstrated through application to several different Ca2+ signaling experiment types.

Simulation of Calcium Dynamics in Realistic Three-Dimensional Domains.

The cytosolic concentration of free calcium ions ([Ca2+]) is an important intracellular messenger in most cell types, and the spatial distribution of [Ca2+] is often critical. In a salivary gland acinar cell, a polarised epithelial cell, whose principal function is to transport water and thus secrete saliva, [Ca2+] controls the secretion of primary saliva, but increases in [Ca2+] are localised to the apical regions of the cell. Hence, any quantitative explanation of how [Ca2+] controls saliva secretion must take into careful account the spatial distribution of the various Ca2+ sources, Ca2+ sinks, and Ca2+-sensitive ion channels. Based on optical slices, we have previously constructed anatomically accurate three-dimensional models of seven salivary gland acinar cells, and thus shown that a model in which Ca2+ responses are confined to the apical regions of the cell is sufficient to provide a quantitative and predictive explanation of primary saliva secretion. However, reconstruction of such anatomically accurate cells is extremely time consuming and inefficient. Here, we present an alternative, mostly automated method of constructing three-dimensional cells that are approximately anatomically accurate and show that the new construction preserves the quantitative accuracy of the model.

A Mathematical Model of Salivary Gland Duct Cells.

Saliva is produced in two stages in the salivary glands: the secretion of primary saliva by the acinus and the modification of saliva composition to final saliva by the intercalated and striated ducts. In order to understand the saliva modification process, we develop a mathematical model for the salivary gland duct. The model utilises the realistic 3D structure of the duct reconstructed from an image stack of gland tissue. Immunostaining results show that TMEM16A and aquaporin are expressed in the intercalated duct cells and that ENaC is not. Based on this, the model predicts that the intercalated duct does not absorb Na[Formula: see text] and Cl[Formula: see text] like the striated duct but secretes a small amount of water instead. The input to the duct model is the time-dependent primary saliva generated by an acinar cell model. Our duct model produces final saliva output that agrees with the experimental measurements at various stimulation levels. It also shows realistic biological features such as duct cell volume, cellular concentrations and membrane potentials. Simplification of the model by omission of all detailed 3D structures of the duct makes a negligible difference to the final saliva output. This shows that saliva production is not sensitive to structural variation of the duct.

Highly localized intracellular Ca2+ signals promote optimal salivary gland fluid secretion.

Salivary fluid secretion involves an intricate choreography of membrane transporters to result in the trans-epithelial movement of NaCl and water into the acinus lumen. Current models are largely based on experimental observations in enzymatically isolated cells where the Ca2+ signal invariably propagates globally and thus appears ideally suited to activate spatially separated Cl and K channels, present on the apical and basolateral plasma membrane, respectively. We monitored Ca2+ signals and salivary secretion in live mice expressing GCamp6F, following stimulation of the nerves innervating the submandibular gland. Consistent with in vitro studies, Ca2+ signals were initiated in the apical endoplasmic reticulum. In marked contrast to in vitro data, highly localized trains of Ca2+ transients that failed to fully propagate from the apical region were observed. Following stimuli optimum for secretion, large apical-basal gradients were elicited. A new mathematical model, incorporating these data was constructed to probe how salivary secretion can be optimally stimulated by apical Ca2+ signals.

Calcium Dynamics and Water Transport in Salivary Acinar Cells.

Saliva is secreted from the acinar cells of the salivary glands, using mechanisms that are similar to other types of water-transporting epithelial cells. Using a combination of theoretical and experimental techniques, over the past 20 years we have continually developed and modified a quantitative model of saliva secretion, and how it is controlled by the dynamics of intracellular calcium. However, over approximately the past 5 years there have been significant developments both in our understanding of the underlying mechanisms and in the way these mechanisms should best be modelled. Here, we review the traditional understanding of how saliva is secreted, and describe how our work has suggested important modifications to this traditional view. We end with a brief description of the most recent data from living animals and discuss how this is now contributing to yet another iteration of model construction and experimental investigation.

A Multicellular Model of Primary Saliva Secretion in the Parotid Gland.

We construct a three-dimensional anatomically accurate multicellular model of a parotid gland acinus to investigate the influence that the topology of its lumen has on primary fluid secretion. Our model consists of seven individual cells, coupled via a common lumen and intercellular signalling. Each cell is equipped with the intracellular calcium ([Formula: see text])-signalling model developed by Pages et al, Bull Math Biol 81: 1394-1426, 2019. https://doi.org/10.1007/s11538-018-00563-z and the secretion model constructed by Vera-Sigüenza et al., Bull Math Biol 81: 699-721, 2019. https://doi.org/10.1007/s11538-018-0534-z. The work presented here is a continuation of these studies. While previous mathematical research has proven invaluable, to the best of our knowledge, a multicellular modelling approach has never been implemented. Studies have hypothesised the need for a multiscale model to understand the primary secretion process, as acinar cells do not operate on an individual basis. Instead, they form racemous clusters that form intricate water and protein delivery networks that join the acini with the gland's ducts-questions regarding the extent to which the acinus topology influences the efficiency of primary fluid secretion to persist. We found that (1) The topology of the acinus has almost no effect on fluid secretion. (2) A multicellular spatial model of secretion is not necessary when modelling fluid flow. Although the inclusion of intercellular signalling introduces vastly more complex dynamics, the total secretory rate remains fundamentally unchanged. (3) To obtain an acinus, or better yet a gland flow rate estimate, one can multiply the output of a well-stirred single-cell model by the total number of cells required.

A Model of [Formula: see text] Dynamics in an Accurate Reconstruction of Parotid Acinar Cells.

We have constructed a spatiotemporal model of [Formula: see text] dynamics in parotid acinar cells, based on new data about the distribution of inositol trisphophate receptors (IPR). The model is solved numerically on a mesh reconstructed from images of a cluster of parotid acinar cells. In contrast to our earlier model (Sneyd et al. in J Theor Biol 419:383-393. https://doi.org/10.1016/j.jtbi.2016.04.030 , 2017b), which cannot generate realistic [Formula: see text] oscillations with the new data on IPR distribution, our new model reproduces the [Formula: see text] dynamics observed in parotid acinar cells. This model is then coupled with a fluid secretion model described in detail in a companion paper: A mathematical model of fluid transport in an accurate reconstruction of a parotid acinar cell (Vera-Sigüenza et al. in Bull Math Biol. https://doi.org/10.1007/s11538-018-0534-z , 2018b). Based on the new measurements of IPR distribution, we show that Class I models (where [Formula: see text] oscillations can occur at constant [[Formula: see text]]) can produce [Formula: see text] oscillations in parotid acinar cells, whereas Class II models (where [[Formula: see text]] needs to oscillate in order to produce [Formula: see text] oscillations) are unlikely to do so. In addition, we demonstrate that coupling fluid flow secretion with the [Formula: see text] signalling model changes the dynamics of the [Formula: see text] oscillations significantly, which indicates that [Formula: see text] dynamics and fluid flow cannot be accurately modelled independently. Further, we determine that an active propagation mechanism based on calcium-induced calcium release channels is needed to propagate the [Formula: see text] wave from the apical region to the basal region of the acinar cell.

A Mathematical Model of Fluid Transport in an Accurate Reconstruction of Parotid Acinar Cells.

Salivary gland acinar cells use the calcium ([Formula: see text]) ion as a signalling messenger to regulate a diverse range of intracellular processes, including the secretion of primary saliva. Although the underlying mechanisms responsible for saliva secretion are reasonably well understood, the precise role played by spatially heterogeneous intracellular [Formula: see text] signalling in these cells remains uncertain. In this study, we use a mathematical model, based on new and unpublished experimental data from parotid acinar cells (measured in excised lobules of mouse parotid gland), to investigate how the structure of the cell and the spatio-temporal properties of [Formula: see text] signalling influence the production of primary saliva. We combine a new [Formula: see text] signalling model [described in detail in a companion paper: Pages et al. in Bull Math Biol 2018, submitted] with an existing secretion model (Vera-Sigüenza et al. in Bull Math Biol 80:255-282, 2018. https://doi.org/10.1007/s11538-017-0370-6 ) and solve the resultant model in an anatomically accurate three-dimensional cell. Our study yields three principal results. Firstly, we show that spatial heterogeneities of [Formula: see text] concentration in either the apical or basal regions of the cell have no significant effect on the rate of primary saliva secretion. Secondly, in agreement with previous work (Palk et al., in J Theor Biol 305:45-53, 2012. https://doi.org/10.1016/j.jtbi.2012.04.009 ) we show that the frequency of [Formula: see text] oscillation has no significant effect on the rate of primary saliva secretion, which is determined almost entirely by the mean (over time) of the apical and basal [Formula: see text]. Thirdly, it is possible to model the rate of primary saliva secretion as a quasi-steady-state function of the cytosolic [Formula: see text] averaged over the entire cell when modelling the flow rate is the only interest, thus ignoring all the dynamic complexity not only of the fluid secretion mechanism but also of the intracellular heterogeneity of [Formula: see text]. Taken together, our results demonstrate that an accurate multiscale model of primary saliva secretion from a single acinar cell can be constructed by ignoring the vast majority of the spatial and temporal complexity of the underlying mechanisms.

Modeling calcium waves in an anatomically accurate three-dimensional parotid acinar cell.

We construct a model of calcium waves in a three-dimensional anatomically accurate parotid acinar cell, constructed from experimental data. Gradients of inositol trisphosphate receptor (IPR) density are imposed, with the IPR density being greater closer to the lumen, which has a branched structure, and inositol trisphosphate (IP3) is produced only at the basal membrane. We show (1) that IP3 equilibrates so quickly across the cell that it can be assumed to be spatially homogeneous; (2) spatial separation of the sites of IP3 action and IP3 production does not preclude the formation of stable oscillatory Ca2+ waves. However, these waves are not waves in the mathematical sense of a traveling wave with fixed profile. They result instead from a time delay between the Ca2+ rise in the apical and basal regions; (3) the ryanodine receptors serve to reinforce the Ca2+ wave, but are not necessary for the wave to exist; (4) a spatially independent model is not sufficient to study saliva secretion, although a one-dimensional model might be sufficient. Our results here form the first stages of the construction of a multiscale and multicellular model of saliva secretion in an entire acinus.